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moremanent magnetization is found to decay algebraically and the temperature de- 
pendent exponents agree very well with the experimentally determined values. The 
nonequilibrium autocorrelation function C(t, t w ) shows a crossover at the waiting (or 
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faster algebraic decay for t^>t w with an exponent similar to one for the remanent 
magnetization. 
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The measurement of dynamical nonequilibrium quantities in real spin glasses has a 
long history. The typical experiment that has been conducted many times is the follow- 
ing: Within a magnetic field the spin glass (for instance Cu(Mn), Au(Fe), Feo.sMno.sTiOs, 
(Fea;Ni(i_ x ))75Pi 6 B 6 Al3, Cd x Mn(i_ x )Te, etc.) is cooled down to temperatures below the 
freezing temperature T g and either immediately or after a certain waiting time t w the field 
is switched off. Then the so called (thermo)remanent magnetization M rem (t) is measured 
as a function of time t. The asymptotic time dependence of this quantity is found to be 
algebraic well below T g (T/T g < 0.98) in the short range Ising spin glass Feo.sMno.sTiOs f||| 
and in an amorphous metallic spin glass (Fe x Ni( 1 _ a ,))75Pi6B6Al3 |J. Furthermore the time- 
dependence of the remanent magnetization depends on the waiting time t w , a phenomenon 
called aging 

Several attempts have been made to explain this behavior theoretically (0,§,^|,1O|,11 and 
a wide variety of functional forms for the time dependence of the remanent magnetization 
is found. The problem lays in the fact that starting from a microscopic model or model- 
Hamiltonian one encounters insurmountable difficulties in trying to solve the nonequilibrium 
dynamics. Therefore additional assumptions have to be made and the final outcome - 



stretched exponential J7|, algebraic [ |TTj] , [TT] ] or logarithmic |5| decay — depends on them. 
Even within the mean field approximation it is hard to obtain any analytical [l2|jll| or 
semianalytical 0] results. 

Once a microscopic model for a spin glass has been formulated, one can in principle try 
to extract its macroscopic behavior via Monte Carlo (MC) simulations. In contrast to an 
analytical treatment, where the calculation of dynamical nonequilibrium quantities within 
the spin glass phase (instead of those characterizing equilibrium, see |15|) is even more com- 
plicated, MC simulation can be done with less effort for certain nonequilibrium situations, 
since equilibration times reach astronomical values in case of spin glasses |T6|Ji7| , p!8[] . The 
remanent magnetization with zero waiting time has been investigated numerically for the 
mean field version of a spin glass model PJT9|| and for the three dimensional EA (Edwards- 
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Anderson) model only right at the critical temperature [20[|. Quite recently attempts have 



been made to investigate the whole temperature range below T g numerically [^T|,^]. How- 
ever, a systematic MC study of nonequilibrium correlations and aging phenomena within 
the frozen phase (T < T c ) of the three-dimensional EA spin glass model has not been made 
up to now. The results such an investigation, its theoretical implications and comparison 
with experiments will be reported in this letter. 

The system under consideration is the three-dimensional Ising spin-glass with nearest 
neighbor interactions and a discrete bond distribution. Its Hamiltonian is 

H = -^JijOiOj , (1) 

where the spins cr, = ±1 occupy the sites of a L x L x L simple cubic lattice with periodic 
boundary conditions and the random nearest neighbor interactions Jy take on the values 
+1 or —1 with probability 1/2. We consider single spin flip dynamics and used a special, 
very fast implementation of the Metropolis algorithm on a Cray YM-P (see ref. p3| for 
details). The simulations were done in the frozen phase, that means at temperatures below 
T c = 1.175 ± 0.025 (see ref. [0,|l]],[il| ) . All measured quantities are averaged over at least 



128 samples (smaller system sizes were averaged over up to 1280 samples). The system size 
was increased until no further size dependence of the results has been observed within the 
simulation time (t < 10 6 ), which is measured in MC sweeps through the whole lattice. It 
turns out that L = 32 is large enough for this time range (cf. |p"8f). 

The system was prepared in a fully magnetized initial configuration and then the sim- 
ulation was run for a time t w (=waiting time) and then the spin configuration g_(t w ) was 
stored. From now on after each MC step (data were then averaged over appropriate time 
intervals, cf. |18|]) the following correlation function was measured: 
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C(t,t w ) = — E(^(* + ^K(^)> > (2) 

i 

where (• • ■) means a thermal average (i.e. an average over different realizations of the ther- 
mal noise, but the same initial configuration) and the bar means an average over different 



realizations of the bond-disorder. The quantity C(t,0) corresponds to the remaining mag- 
netization of the system after a time t. 

M taa (t) = C(t,0) . (3) 

This quantity is directly related to the experimentally determined thermoremanent magneti- 
zation with zero waiting time (i.e. without aging) and nearly saturated initial magnetization. 
The result for the remanent magnetization M rem (t) is shown in fig. |l] within a log-log plot. Its 
decay clearly obeys a power law for large times and temperatures in the range 1.1 > T > 0.5. 
The exponent A(T) for the fit 

M rem (t) ex r x ^ , (t > 10 2 ) (4) 

is plotted in fig. 0, upper curve. It starts at A = 0.36 ± 0.01 for T = 1.1 (and can be 
extrapolated via the fit indicated in figure 2 to 0.39±0.01 for T=T C , which was already 
found in ||20|| ) and decreases monotonically with temperature. For the short range Ising 
spin glass Feo.sMno.sTiOs and for certain amorphous metallic spin glasses not only the same 
algebraic decay of the remanent magnetization has been observed, but also the shape of the 
functional temperature dependence of the exponent A(T) and even its numerical values are 
in excellent agreement: from figure 4b in ref. f| one may for instance read off A(T 9 ) « 0.38 
and A(0.5T g ) « 0.12, concurring within the errorbars to the corresponding data plotted in 
fig. 2. 

It was argued || that the decay of M rem (t) should be logarithmic (i.e. M Tem (t) oc 
(lnt) _A// ^) below T c , but a fit of the data in fig. 1 for T > 0.5 does not yield accept- 
able results over the range of the observation time. We want to focus some attention to the 
T = 0.4 curve: It bends upward in the log-log plot for t > 10 4 , which could indicate the 
onset of a slowlier than algebraic decay for temperatures smaller than 0.5, logarithmic for 
instance. 

Another indication that something new happens at lower temperatures can be obtained 
by looking at the short time behavior of M rem (t): At T « 0.5 a plateau begins to develop 
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for t < 10 2 , which can clearly be seen for T = 0.4 and becomes even more pronounced and 
wider for even smaller temperatures. It can be excactly reproduced in shape and location 
for smaller and larger sizes and number of samples, which means that it is a physical effect 
and not only a fluctuation. A possible interpretation might be that the system gets trapped 
in metastable states, whose lifetime grows with decreasing temperatures. 

Next we turn our attention to the correlation function C(t,t w ) with t w = 10 a (a = 
1, . . . , 5). In contrast to M rem (i) the correlation function C(t, t w ) for t w ^ is not directly re- 
lated to the thermoremanent magnetization M(t, t w ) at time t + t w in a temperature-quench 
experiment, where the field H is switched off at time t w after the quench (see e.g. ||24j| ). In 
equilibrium (t w — > oo) C(t, oo) is related to the relaxation function R(t, oo) = M(t,oo)/H 
via the fluctuation dissipation theorem (FDT) R(t, oo) = C(t, oo)//c B T. However, in a 
nonequilibrium situation, like the one considered here, slight differences between them exist 
pl| , |25[| (e.g. in the location of the maximum relaxation rate). The magnetization that is 
induced by a small external field (H <C 1) for model (|1|) is rather small, therefore the func- 
tional form of M(t,t w ) is harder to determine accurately via MC-simulations. This is the 
reason why in this letter the focus is on C(t,t w ). 

A typical set of data for a particular temperature (T = 0.8) is shown in fig. ^ in a log-log 
plot. One observes a crossover from a slow algebraic decay for t <C t w to a faster algebraic 
decay for t ^> t w . The crossover time is simply defined as the intersection of the two straight 
line fits for short- and long-time behavior in the log-log plot. For the long-time behavior 
the fit to 

C{t, t w ) oc r X( - T ' tw) , t^>t w (5) 

yields a set of exponents that is depicted in fig. ^. For increasing t w the exponent X(T,t w ) 
decreases only slightly and the waiting time dependence becomes weaker for lower temper- 
atures. By looking at fig. one observes that it is difficult to extract X(T,t w ) for t w = 10 4 
and 10 5 since there are only 2 or 1 decades left to fit the exponent — therefore they are not 
shown in fig. 2. The exponent describing the short time (t t w ) behavior of C(t,t w ), 
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C(t,t w ) oc/T^ , t<t w , (6) 

which is depicted in fig. f|, is independent of the waiting time t w . Since the system was able 
to equilibrate over a time t w , all processes occuring on timescales smaller than t w have the 
charcteristics of equilibrium dynamics and therefore the exponent x(T) is identical to that 
describing the decay of the equilibrium autocorrelation function q(t) = lim^-nx, C(t, t w ). 
The latter was investigated in |l8j and the exponents that are reported there for T > 0.7T C 
agree with the values shown in fig. ||. They also agree with those determined experimentally 
in the short range Ising spin glass Feo.sMno.sTiOs via the above mentioned relaxation 
function R(t, t w ) = M(t, t w )/H for t <^t w (note that in this quasiequilibrium-regime C (t, t w ) 
and R(t,t w ) are related via the FDT |2l|j25|] , yielding the same exponents for both): close 
to T g {T/Tg = 1.029) they obtain x = 0.07. Furthermore there seems to be a temperature 
at about 0.3, where x(T) becomes zero, which could be another indication for the above 
mentioned onset of a logarithmic decay of the correlation functions . 



Although the decay of the nonequilibrium correlations in the temperature range 0.5 < 
T < 1.1 is algebraic rather than logarithmic as predicted by the droplet picture proposed in 
||, this picture might not be inappropriate: Let us assume the following scaling law for the 
dependence of the free energy barriers B on a length scale L of the regions to be relaxed: 
B oc AlnL instead of B oc as in fl9|. Then one ends up with an algebraic decay of e.g. 
the remanent magnetization by observing (see 0) that the typical length scale of domains 
R t now grows with time according to A\nR t ~ Tint, which means R t oc t T / A ( T ), leading to 
equations (0-^1). 

In the context of the phenomenological model for the dynamics and aging in disordered 
systems developped in ref . |LT| , the algebraic decay of correlations found so far implies that 



the probability distribution of free energy barriers is exponential in the temperature range 
of 0.5 < T < 1.1 for the system under consideration. Furthermore we would like to mention 
that a fit to the functional form for the short time behavior (t <ti t w ) C(t, t w ) ~ 1 — a(t/t w ) y 



proposed in [0 works also quite well for our data, although not as convincingly as equation 
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Guided by equations and (Q) we tried to put our results into the following scaling 
form: 

C(t,t w ) = err*™ $ T {t/t w ) , (7) 

where <&t(u) = 1 for y = and $r(y) oc y 1 ^"^^ for y — > oo. The form (0) has re- 
cently been used |27j successfully to extract the critical dynamical exponent z from the 



nonequilibrium correlation function (0) via finite size scaling, where the waiting time t w 
has been replaced by the relaxation time r oc L z in the critical region. For temperatures 
below T = 0.8 equation (|7|) yields an acceptable fit (which can already be deduced from the 
neglegible waiting time dependence of the exponents \(T,t w ) for T < 0.7, see fig. 2). 

Concluding we have reported new results of numerical nonequilibrium simulations 
that show an excellent concurrence with experiments on the short range Ising spin glass 
Feo.5Mno.5TiO"3 and on amorphous metallic spin glasses: not only a single exponent but a 
whole continuum of (temperature dependent) exponents for the remanent magnetization are 
found to agree within the numerical errors. Although the values for the exponents extracted 
from experiments might vary somewhat depending on the microscopic details (range of inter- 
actions, spin-type) the main features of the relaxation and the dynamics of many different 
three-dimensional spin glasses are very similar and the functional forms of the remanent 



magnetization decay should be the same for different systems [^8 . 

Furthermore we have shown that aging phenomena in the spin glass model under con- 
sideration can be observed via the measurement of a particular correlation function and 
that its nonequilibrium dynamics is indeed gouverned by its equilibrium characteristics for 
time scales smaller than the imposed waiting (or aging) time. This gives an interesting new 



perspective (see also [^,^]) to extract equilibrium quantities, which are hard to obtain via 
MC simulations within the spin glass phase. Finally, by observing plateaus in the short 
time behavior and slowing down of the algebraic decay of the remanent magnetization, we 
revealed a dynamical scenario at very low temperatures that is not yet fully understood. 
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FIGURES 



FIG. 1. The remanent Magnetization M rem (i) versus the time i in a log-log plot for varying 
temperatures. From top to bottom we have T=0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 and 1.1. The sytem 
size is L=32 and the data are averaged over 128 samples. The errorbars are of the size of the 
symbols for M rem < 0.01 and much smaller for larger M rem . 

FIG. 2. Upper curve: The exponent A(T) for the remanent magnetization (||) versus tem- 
perature. The points represented by diamonds (o, plus errorbars) are those extracted from figure 
| and the full curve is a least square fit to a quadratic polynomial as a guideline to the eye. 
Lower Points: The nonequilibrium exponent X(T,t w ) (see equation (||)) extracted from the 
long-time behavior (t S> t w ) of the nonequilibrium correlation function C(t,t w ) (see figure ||) 
for fixed values of t w versus the temperature T. From top to bottom we have: (A) t w =10, (□) 
t w =100 and (o) £^=1000. The errorbars are indicated. 

FIG. 3. The averaged nonequilibrium spin autocorrelation function C(t,t w ) of equation (||) for 
fixed values of t w versus time t on a double logarithmic time-scale. The temperature is fixed to be 
T=0.8 and from top to bottom we have i^=10 5 , 10 4 , 10 3 , 10 2 and 10. The system size is L=32 
and the data are averaged over 128 samples. The size of the errorbars is only a fraction of the size 
of the symbols. 

FIG. 4. The equilibrium exponent x(T) for the equilibrium autocorrelation funcion q(t) ex- 
tracted from the short-time behavior (t <C t w ) of the nonequilibrium correlation function C(t,t w ) 
(see equation (||)) versus the temperature T. The errorbars are smaller than the circles, as in- 
dicated. For T < 0.8 the data are fitted to a straight line, which shows that at approximately 
T = 0.3 the exponent x(T) vanishes. 
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